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Abstract 

We study the correlations between eigenvalues of the large random matri- 
ces by a renormalization group approach. The results strongly support the 
universality of the correlations proposed by Brezin and Zee. Then we apply 
the results to the ground state of the ID particles with long range interac- 
tions in a confinement potential. We obtain the exact ground state. We also 
show the existence of a transition similar to a phase separation. Before and 
after the transition, we obtain the density-density correlation explicitly. The 
correlation shows nontrivial universal behavior. 
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The random matrix theory has a long and distinguished history |l|] and the results 
have been applied to many fields, including heavy nuclei, quantum gravity and mesoscopic 
systems [|. 

Recently Brezin and Zee 0] ||^ have proposed the universal correlations between eigen- 
values for a broad class of large random matrices by calculating them explicitly with an 
ansatz for the orthogonal polynomials. What is remarkable in their results is that, except 
for an overall scale setting, the width of the correlations is completely independent of the de- 
tails of the Hamiltonians. Beenakker ^ extended their results by a different technique. As 
is suggested by Brezin and Zee, we try to understand this through a renormalization group 
(RG) approach. Following the RG method developed by Brezin and Zinn- Justin 0, we 
calculate the /3 function for a broad class of unitary ensembles and show the existence of the 
stable Gaussian fixed point. This result strongly supports the existence of the universality 
in the large random matrices proposed by Brezin and Zee. 

Moreover we consider the ID particles with long range interactions in a confinement 
potential using the random matrix theory. The exact ground state is obtained. We also 
obtain the density distribution and the density-density correlation explicitly and investigate 
them. The ground state shows a novel transition to the "separation" phase as described 
below. 

Let us consider the random matrix model described by the Hamiltonian 

H=\TrM^+ E [^(TrM^)^ + ^(TrM^^)], (1) 

where M is an by Hermitan matrix. Set one of pi^'s or qk^s to be 5f( 7^ ) and the rest 
to be zero. The free energy of this model is F{g) = — limAr^oo dM exp{—H)], where 

dM is AfJJi dMii Y[i<j dReMijdlmMij and Af is set to satisfy F{g = 0) = |. Then we have 
the differential equation Q [|^ 

CO o 771 00 ^2 p 00 

(E «^^'')(^) = -^(E bK9''){%;^) + (E CKg""), (2) 

K=0 K=0 K=0 

and the P function is given by [|| 
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^^^^ = a^q^ 

Equation (2) with (3) denotes the large- hmit of the Callan-Symanzik-hke equation. The 
idea is an analogy between the matrix size in the matrix model theory and the momentum 
cut off in the quantum field theory, which the double scaling limit in the 2D quantum gravity 
suggests. The /3 function has two zero points, one is the Gaussian fixed point {g = 0) and 
the other is the fixed point {g < 0) which governs the double scaling limit in the 2D quantum 
gravity 0. The latter fixed point, however, does not play any role in the present problem 
since the coupling constant g is negative. 

The crucial step for understanding the universal correlations is to investigate the stability 
of the Gaussian fixed point. In the following we put ao = 1 without loss of generality. Hikami 
investigated the case when q2 is g and obtained 

Pi9) = -4i±i^ (4) 
^^^^ 3(1 + 24^). ^ ^ 

This is consistent with the known exact results obtained by different techniques 0. For 

example, (4) implies that a phase transition takes place ai g = while the exact critical 

value of g is also — ^. 

If we have the complete series expansion of F{g), it is possible to calculate the f3 function 
in principle. In general, however, it is technically difficult to obtain high order terms. Thus 
we shall calculate the series expansion up to low order and determine the values of ai , bo 
and bi. We have set ax = bx = ck-i = for K = 2, 3, ■ ■ ■ . They are still expected to be 
good approximations, as demonstrated in several cases §]. 

We consider two classes of models. 

i) i/= iTrM2 + ^(TrM2f. 

At first we calculate the series expansions of F{g). F{g) is 

1 f n1 1 g 

^ = ~ i'^o.N^^''^J '^^''^ " e^P(-2^-]^^")]' (5) 

where x=TrM'^ and c is a constant. By replacing x by N'^x, we obtain the free energy by the 
saddle point method as F = jx+gx"^ — ^ \nx. The saddle point equation is ^x+ngx'"- — j = 0. 
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Thus we have |^ = l-2n'^g+0{g'^) in the small g limit and |^ = ^n+i n+i +0{-^) 

in the large g limit. 

Putting these series expansions into (2), we obtain the (3 function from (3) as 

^ g{l + 2n^g) 

n[l + {2n^ + 2n)g]. 

When n=2, this is the exact solution of (2) and (3) and consistent with the exact results 
obtained by other techniques |T^. Thus we expect that the (3 function for n=2 is exact. 



ii) H =1 Tr + ^ Tr M^". 

Since the matrix M can be diagonalized by a unitary matrix, the free energy can be 
expressed by the eigenvalues. Thus one can have the complete series expansion of F{g) 
for n=2 and the /3 function is obtained as (4). It is, however, difficult to get the 
exact results for higher n. Instead we exactly map the Hamiltonian to H = ^ Tr + 
J2k=i jv2nK-2 (Tr M^)"-'^, integrating out the angular part of the Hamiltonian [Q. In the 
lowest order in g, the Hamiltonian is approximated by = | Tr + jy2n-2 (Tr M^)", where 
di is obtained by comparing F before and after the mapping. The result is di = 2{n = 
2),5(n = 3), lA{n = 4),42(n = 5), 132(n = 6), ■ ■ ■. By following the same procedure as (i), 
we obtain 

Pig) = 9il + 2n^d,g) 

n[l + {2n'' + 2n)dig]. ^' 

The results (6) and (7) show that the Gaussian fixed point (^'=0 ) is stable for a broad 

class of large random matrices ( see Fig. 1 ) . This simple renormalization group argument 

leads to a direct and intuitive understanding of the universal properties of the large random 

matrices . 

Next we shall apply the results obtained above to ID A^-body Fermi or Bose systems in 
a confinement potential. The Hamiltonian is 

H = Hj + Hn + Hui, (8) 
= -E3Z2 +E7— (9) 



and 
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(10) 



Hjn = Y.Dix,-x^f, (11) 

where Hj is the Hamiltonian of the A^-body system with ^ interaction, Hjj is a confinement 
potential and Hjjj is attractive harmonic interaction. Note that Hu is double-well like when 
B is negative and C is positive. Let us parameterize A, B, C, D and qq by 

A = - 12« - 12(iV - 1)«A, (12) 

B = Sau, (13) 

C = lQa\ (14) 

D = 4aX, (15) 

go = 2X' - 2A, (16) 

where we set a to be positive. Although the parameters A, B, C, D and qq are not totally 
independent, we can change the depth of the double well without changing the interaction 
by changing u;, while keeping a and A fixed, where u can also be negative and the new uni- 
versality class described below appears. At a = , we have the Calogero-Moser-Sutherland 



(CMS) model |jT2[ |[T^ , which is known to be integrable and has been studied ex- 

tensively. The system described hj H = Hj + Hjjj was also described by Calogero [jl2|. 
Thus we call (8) with (9), (10) and (11) the generalized Calogero-Moser-Sutherland (GCMS) 
model. We can also generalize Hu and Hm in a way that the extended Hamiltonian still 
has the universal properties investigated below, but we shall restrict our attention to (8) for 
simplicity. 

It can be shown that the following Jastrow-type function is an exact eigenfunction with 

E = uN[l + X{N -1)] 

I N N 

"^{xi, ■■■,xn) = SY[\xi- Xj|^exp(--u;^x- - a^x-). (17) 

i<j i=l i=l 



The factor S is equal to 1 for bosons, while for fermions S = (—1)'^ is the parity of the 
particle ordering permutation P, xp^^ < xp^ < ■ ■ ■ < xp^. It is the ground state since, for 



each ordering of the particles, this wavefunction is nodeless [|Tl| |[T^. Thus we obtain the 
exact ground state. Note that even for the much simpler system of non-interacting particles 
in an anharmonic potential, the ground state can not be obtained exactly. The squared 
amplitude of the wavefunction is 

N N 

■ ■ ■ ,X7v)p = J]^ |xi - Xjf exp{-ujJ2^1 - 2a Xl^^)' W 

i<j 1=1 1=1 

where (3 is 2 A. 

For /3=1, 2 and 4, (18) is identical to the probability weight of the orthogonal, unitary 
and symplectic ensembles in the random matrix theory respectively, where we identify the 
eigenvalues of random matrices with the particle coordinates. Thus the following arguments 
are also applicable to the random matrix theory. 

The exact ground state (18) has nontrivial structures. In order to investigate them we 
take the limit — oo, — ~ 0(1) and ^ ~ taking into account the 

effect of the confinement potential. At first, we consider the density distribution p{x). In 
the large- limit we can apply the saddle point method to (18) and obtain 



uox + Aaox' = ^P I dy^^ (19) 
2^ J X — y, 



where P denotes the principal part of the integral. This integral equation can be sloved by 
the Hilbert transformation 0] and we have (see Fig. 2), 
i) for luq > ujc 

\ ^(^^2 + ^ + ^){c? - x2)i (-a < X < a) 



(otherwise) 



ii) for uq < LOc 



p{x) = <) - /3 ' ^ (21) 

(otherwise) , 
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with a- = ^[-^ + + i^)^],c^ = -^- ,f^,d^ = -il + x/£ and 



24ao I- /3 ^ p J 11 4ao y 4ao ' 4ao V 4ao 

-74^. 

When uo < Uc, the paticles become separated (see Fig. 2). At u>c, the universahty class 
of the system changes. We call this novel phenomenon "separation" . A similar phenomenon 
is the phase separation in the t-J model |14]. 



Secondly, we consider the density-density correlation of the GCMS model. As a direct 
consequence of the Brezin and Zee's results and the extended results by Beenakker 

0, we have 

(23) 

where \x — y\ ~ 0(1) and Pc{x, y) is defined by < J2i — Xi) J2j ^iu — Xj) > — <J2i ^x — 
Xi) >< J2j S{y — Xj) >. Notice that there is the effect of the confinement potential in (22), 
since (22) is not the function of only {x — y). The importance of this result is that Pc{x, y) is 
universal in contrast to p{x) . Although (22) is correct for ujq > Uc, the numerical simulation 
|T3| suggests that it does not hold for ujq < tOc, where the potential has a deep double well. 
We expect that the system belongs to a different universality class for ujq < Uc- In order 
to show this explicitly, we calculate Pc{x, y) after separation using the method developed 
recently by Beenakker 0. Let us consider the following wavefunction in which the support 
of p{x) is (p, g) U (r, s) {p < q < r < s) (see Fig. 3) 

\^{xi,--- ,xn)\^ = exp[/3^1n|x, -XjI -/3^y(xi)], (24) 

i<j i 

where V{x) is a function. In the large-iV limit we can apply the saddle point method to 
(24) and obtain 



P 



rdyJ^ + PrdyJ^='^. (25) 
Jp X — y Jr X — y dx 



A variation of (25) gives P Jp dy^^^ + P dy^^^ = Since separation does take place 

and there are no particles in the finite range between q and r, we have 



5 I dxp{x) = S f dxp{x) = 0. (26) 

Thus, from the fact that pc{x,y) is symmetric under an exchange of x and y, and the 
Beenakker formula pc{x,y) = —^jp^ we have 



Pc{x,y) 



(otherwise) , 

where |a; — y| ~ 0(1) and Aq and So are constants. Note that Pc{x, y) is still universal after 
separation. 

In summary, we calculated the (3 functions for a broad class of random matrices and 
found that the Gaussian fixed point is stable . As is suggested by Brezin and Zee @] 
the stable Gaussian fixed point is responsible for the universality of the correlations 
in the large random matrices. This result is interpretated in the GCMS model (8). We 
investigated the system in the large limit. The system can be in the two universality classes 
depending on the parameters of the Hamiltonian. We call this change of the universality class 
"separation" . This phenomenon is similar to the phase separation in the t-J model. In all the 
parameter regions , we obtained the density distribution and the density-density correlation 
explicitly. The density-density correlation shows the universal behavior in contrast to the 
density distribution. We take into account the effect of the confinement potential and our 
results may have some implications to the mesoscopic systems such as quantum wire. 

We wish to thank S. Hikami, T. S. Kobayakawa, and A. Zee for very helpful discussions. 
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FIGURES 



Fig. 1. The (3 function of the random matrix model (schematic). 

The arrow denotes the renormahzation group flow . There are two fixed points. One is 
the stable Gaussian fixed point and the other is the unstable fixed point which governs the 
double scaling limit in the 2D gravity 0] . 

Fig. 2. The density distribution of the generalized Calogero-Moser-Sutherland model 
(schematic). 

X denotes the space coordinate and p{x) the density distribution. 
Fig. 3. The density distribution of the wavefunction (24) (schematic). 
X denotes the space coordinate and p{x) the density distribution. 



10 



Fig.2 
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